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Abstract 

In this work we study the process e+ + e~ — > J/^ + r]^ at energy ^/s = 10.6Gey observed 
recently at B-factories whose measurements were made by Babar and Belle groups. We calculate 
the cross section for this process in the Bethe-salpeter formalism under Covariant Instantaneous 
Anstaz (CIA). To simplify our calculation, the heavy quark approximation is employed in the quark 
and gluon propagators. In the exclusive process of e"'"e~ annihilation into two heavy quarkonia, 
the cross section calculated in this scenario is compatible with the experimental data of Babar and 
Belle. 
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I. INTRODUCTION 

In this work we study the exclusive production process — > J/tp+t]^ at energy ^/s = 
lO.GGeV observed at B-factories [1,2,3] whose measurements have recently been done by 
Babar and Belle groups. It is well known that there was a significant discrepancy between the 
experimental measurements P-|3] and the non-relativistic QCD (NRQCD)[11 15] predictions 
for this process at centre of mass energies ^/s ~ lO.GGeV. This process has been recently 
studied in a Bethe-Salpeter formalism [6] in Instantaneous Approximation (lA). To simplify 
calculations, the authors have employed heavy quark limit in the propagators for studying 
systems composed of heavy charm and anti-charm quarks. 

We wish to mention that Bethe-Salpeter equation (BSE) is a conventional non- 
perturbative approach in dealing with relativistic bound state problems in QCD. It is firmly 
established in the framework of field theory and from the solutions we can obtain useful 
information about the inner structure of hadrons, which is also crucial in treating high en- 
ergy hadronic scattering and production processes. Despite its drawback of having to input 
model-dependent kernel, these studies have become an interesting topic in recent years, 
since calculations have shown that BSE framework using phenomenological potentials can 
give satisfactory results as more and more data is being accumulated. Further by adopting 
this framework we get more insight about the treatment of this process. This is mainly 
due to the unambiguous definition of BS wave function which is expressible by time ordered 
product of Heisenberg picture operators. This provides exact effective coupling vertex for 
bound state particle with all its N (N = 2 for mesons) constituents and can be considered 
as summing up all the non-perturbative QCD effects in the bound state. 

On lines of [6], we try to study this process in the framework of BSE under Covariant 
Instantaneous Ansatz (CIA) which is a Lorentz- invariant generalization of Instantaneous 
Approximation (lA). What distinguishes CIA from the other 3D reductions of BSE is its 
capacity for a two-way connection: an exact 3D BSE reduction for a qq system, and an 
equally exact reconstruction of original 4D BSE, the former to make contact with the mass 
spectrum, while the latter for calculation of transition amplitudes as 4D quark-loop integrals 
[7I4TT]. We wish to emphasise here that in these studies one of the main ingredients is the 
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Dirac structure of the Bethe-Salpeter wave function (BSW). The copious Dirac structure 
of BSW was aheady studied by LlewUyn Smith [12j much earher. Recent studies [131 [H] 
have revealed that various covariant structures in BSWs of various hadrons are necessary to 
obtain quantitatively accurate observables. It has been further noticed that all covariants 
do not contribute equally for calculation of meson observables. To address this problem, re- 
cently we thought of investigating how to arrange these covariants systematically in BSWs. 
Thus, in a recent work [lOj, we developed a power counting rule for incorporating various 
Dirac structures in BSW, order-by-order in powers of inverse of meson mass. According to 
this power counting rule,the Leading order (LO) covariants are expected to contribute max- 
imum to calculation of any meson observable, followed by the next-to-leading order (NLO) 
covariants. Taking in view of this fact, as a first step we have outlined the Dirac covari- 
ants and expanded the coefficients to the leading order (LO), and calculated the leptonic 
decay constants of vector mesons [10] as well as pseudoscalar mesons [llj at this order. The 
results were found to be close to data. In another recent work [T3], we studied leptonic 
decay constants of unequal mass pseudoscalar mesons like n,K,D, D^and B and radiative 
decays of equal mass pseudoscalar mesons like 7r°, t]^ by taking into account both the leading 
order (LO) and Next-to-Leading Order (NLO) Dirac covariants. It was found that the con- 
tribution of leading order (LO) covariants to decay constants was maximum (about 90-95 
percent) for heavier mesons composed of c and b quarks like D,Ds,B and 1]^ [15j , while 
there was little contribution from NLO covariants. We now calculate the cross section for 
the process e~ + e'^ ^ J/ip + rj^ by employing the most leading of the LO covariants (such 
as 75 for heavy pseudoscalar mesons like t]^ and ij.e for heavy vector meson like J/ip com- 
prising of heavy charm and anti-charm quarks for which BS formalism is quite suitable. 
In order to simplify the calculation, we will further impose the heavy quark approximation 
(P ~ M, q << M) on the quark and gluon propagators to simplify the integrals as in Ref.[6]. 
Under this approximation, our results are comparable with the data [THSj. The remainder 
of this paper is organized as follows: In sec. II, we will study the BS equations for vector and 
pseudoscalar quarkonia. In sec. Ill, we will calculate the amplitude and cross section for the 
process + e~ — > J/"^ + rj^ in the BS formalism. Sec. IV is reserved for conclusions and 
discussions. 
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II. THE BETHE-SALPETER WAVE FUNCTION UNDER CIA 



We briefly outline the BSE framework under CIA. For simplicity, lets consider a qq system 
comprising of scalar quarks with an effective kernel 4D wave function $(P, g), and with 
the 4D BSE, 

z(27r)^AiA2$(P, q) = j dU'K{q, g'), (1) 

where Ai 2 = ml 2 + Pi 2^^^ inverse propagators, and mi^2 are (effective) constituent 
masses of quarks. The 4-momenta of the quark and anti-quark, pi^2, are related to the 
internal 4-momentum and total momentum of hadron of mass M as Pi,2^ = '"^i,2PAt±g/i, 
where mi,2 = [1 i (^1 ~ m2)/M'^]/2 are the Wightman-Garding (WG) definitions of masses 
of individual quarks. Now it is convenient to express the internal momentum of the hadron 
q as the sum of two parts, the transverse component, q^ = q^ — p?P^ which is orthogonal 
to total hadron momentum P (ie. q.P = regardless of whether the individual quarks are 
on-shell or off-shell), and the longitudinal component, crP^ = {q ■ P/P^)P^, which is parallel 
to P. To obtain Hadron-quark vertex, use an Ansatz on the BS kernel K in Eq. ([!]) which 
is assumed to depend on the 3D variables q^, [3 [HI [10], [Til [IS] 

K{q,q') = K{q,q'), (2) 

(A similar form of the BS kernel was also earlier suggested in ref. [E]). Hence, the 
longitudinal component, aP^ of q^, does not appear in the form K{q, q') of the kernel. For 
reducing Eq.(l) to the 3D form of BSE, we define a 3D wave function 0(g)as, 

/+00 
Mda<^{P,q) (3) 

Substituting Eq.(3) in Eq.(l), with the definition of the kernel in Eq.(2), we get a covariant 
version of the Salpeter equation which is in fact a 3D BSE: 

(27r)3D(g)0(g) = J d^'K{q,q')<P{q'). (4) 

Here D{q) is the 3D denominator function defined below whose value is obtained by 
evaluating contour integration over inverse quark propagators in the complex cx-plane by 
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noting their corresponding pole positions [TOl [HI [T7]as. 

^1,2 = "^1,2 + (5) 



+ 00 1,1 



r Mda _ + 
D{q) 27r2 7 A1A2 + c^a)' - ' "^^'^ 



We note that the RHS of Eq.(4) is exactly identical to the RHS of Eq.(l) by virtue of Eq.(2) 
and Eq.(3). We thus have an exact interconnection between 3D BSE and the 4D BSE, and 
hence between the 3D wave function0(g) and the 4D wave function^ (P. q) [7-nll\ : 

A,A,t(P.,) = ^MM«,rM). (6) 

where T{q) is the Bethe-Salpeter Hadron-quark vertex function for a meson comprising of 
scalar quarks. 

The 4D BS wave function $(P, g) can be reconstructed from the 3D BS wave function 
0(g) as: 

HP,<l) = ^m^, (7) 

where A, = (mf + p^),(i=l,2) are the inverse propagators for scalar quarks which flank 
the hadron-quark vertexF. This 4D hadron-quark vertex r(g) satisfies a 4D BSE with 
a natural off-shell extension over the entire 4D space (due to the positive definiteness of 
the quantity (f' = q^ — {q ■ PY / P"^ throughout the entire 4D space) and thus provides a 
fully Lorentz-invariant basis for evaluation of various transition amplitudes through various 
quark loop diagrams. Due to these properties, this framework of BSE under CIA can 
be profitably employed not only for low energy studies but also for evaluation of various 
transition amplitudes at quark level all the way from low energies to high energies. However 
this 4D hadron-quark vertex F is still unnormalized and can be normalized as will be shown 
in the realistic case of fermionic quarks next. 

For fermionic quarks the BSE under CIA can be written as: 

i{2Tit^{P, q) = SMPi)Sf2{P2) J d'q'K{q, q')^{P, g'), (8) 

where the scalar propagators A^^ in the above equations are replaced by fermionic prop- 
agators Sp- Further the Hqq vertex would be a 4 x 4 matrix in the spinor space for which 
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we should incorporate the relevant Dirac structures. For incorporation of the relevant Dirac 
structures in T{q), they are incorporated order-by-order in powers of inverse of meson mass 
M, in accordance with the power counting rule we developed in jTUj. Our aim of developing 
the power counting rule was to find a "criterion" so as to systematically choose among various 
Dirac covariants from their complete set to write wave functions for different mesons (vec- 
tor mesons, pseudoscalar mesons etc.) p^ ITT] . In another recent work, we studied leptonic 
decay constants of unequal mass pseudoscalar mesons like tt, K, D, DgandB and radiative 
decays of equal mass pseudoscalar mesons like 7r°,r7^ by taking into account both the lead- 
ing order (LO) and Next-to-Leading Order (NLO) Dirac covariants. It was found that the 
contribution of leading order (LO) covariants to decay constants was maximum (about 90- 
95 percent) for heavier mesons composed of c and b quarks like D ,Da,Baii(i r]^^^, while 
there was little contribution from NLO covariants. Among the LO covariants, it was also 
noticed that for pseudoscalar mesons, the contribution from covariant, 75 was maximum and 
similarly for vector mesons, the contribution of LO covariant Z7.£:was maximum. 

Thus to simplify the calculations, as a first step, we calculate the cross section for the 
process + J li' + Vc^Y employing the most leading of the LO covariants such as 

75for heavy pseudoscalar mesons like r^^and i'j.efoT heavy vector meson like J/^/'comprising 
of heavy charm and anti-charm quarks for which BS formalism is quite suitable. The full 
fledged normalized 4D BS wave functions for a qq meson with quarks total momentum Pand 
relative momentum q and with individual quarks with momenta pi and p2 can be written 
as: 



where the 4D BS hadron-quark vertex function which absorbs the 4D BS normalizer N is. 



where Fj = 75,z7.e, ... are the relevant Dirac structures for pseudoscalar mesons, vector 
mesons etc. The 4D BS normalizer which is determined from standard current conserving 
conditions and is worked out in the framework of Covariant Instantaneous Ansatz (CIA) 
to give explicit covariance to the full fledged 4D BS wave function, \I'(P, g) and hence to 
the Hadron-quark vertex function, F(g) employed for calculation of transition amplitudes at 



^{P,q) = Sp{p^)Tiq)SFi-p2) 



(9) 



T{q) = NT,D{q)4>{q)/2m 



(10) 
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high energies. 

Thus the 4D hadron-quark vertex function for rj^ is, 



T^{q,) = j,NpD{q,)(P{q,)/2m, (11) 
while the 4D hadron-quark vertex function for J/^ is, 

r^(g;) = t^.eNvD{qa)(l){qa)/2m. (12) 

Here Np and Ny are the 4D BS normahzers for rj^ and J/\l/ meson (with internal momenta 
qa and qb respectively) which are determined through the current conservation condition 
O m Uni EH [15] , whUe D{qa) and D{qb) are the respective denominator functions. (t>{qa) and 
(p{qb) are the 3D BS wave functions for r]^ and J/ip respectively, while e is the polarization 
vector for the J/tp meson. 

As far as the input kernel K{q, g')[71 El [iDl [HI [E], in BSE is concerned, it is taken as 
one-gluon-exchange like as reg ards color [(A(^V2) ■ (A^^Vs)] and spin {'j^f^'^'y^^^) dependence. 
The scalar function V{q — q') is a sum of one-gluon exchange Vqge and a confining term 
Vconf- Thus we can write the interaction kernel as [71 [TO] : 

Kiq,q') = Qa«) ■ Qa(^)) i;(^)v;(^V(g-g'); 



127r ^ 



r2(l + 4aomim2MV)-i/2_ C| 



g«('J-9')-»'- 



(13) 



The Ansatz employed for the spring constant cj^, in Eq. ( 13 ) is IH [H] [TOj [TTl [15 



u;2 ^ Amim2M>ujlas{Ml), M> = Max{M, mi + ma) (14) 
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Here mi, m2 are the Wightman-Garding definitions of masses of constituent quarks de- 
fined earlier. Here tfie proportionality of on asiQ"^) is needed to provide a more direct 
QCD motivation to confinement. This assumption further facilitates a flavour variation in 



u^g. And lOq in Eq. (13) and Eq. (14) is postulated as a universal spring constant which 
is common to all flavours. Here in the expression for V{q — q'), as far as the integrand of 
the confining term Vconf is concerned, the constant term Cq/uI is designed to take account 
of the correct zero point energies, while Oq term (oq ^ 1) simulates an effect of an almost 
linear confinement for heavy quark sectors (large mi, 7712), while retaining the harmonic 
form for light quark sectors (small mi, 7712) |7J as is believed to be true for QCD. Hence the 
term r^(l + 4aomim2^>'"^)^^^^ in the above expression is responsible for effecting a smooth 
transition from harmonic (qq) to linear (QQ) confinement . The values of basic constants 
are: Cq = .29, coq = 0.158 GeV, m„,d = 0.265 GeV, m, = 0.415 GeV, m^ = 1.530 GeV and 
7nh = 4.900 GeV which have been earlier fit to the mass spectrum of qq mesons [7] obtained 
by solving the 3D BSE under Null-Plane Ansatz (NPA). However due to the fact that the 
3D BSE under CIA has a structure which is formally equivalent to the 3D BSE under NPA, 
near the surface P.q = 0, the qq mass spectral results in CIA formalism are exactly the 
same as the corresponding results under NPA formalism[71 [HI [17]. The details of BS model 
under CIA in respect of spectroscopy are thus directly taken over from NPA formalism (For 
details, see [71 [H HI]). As far as the 3D wave function </)(g) is concerned, it satisfies the 3D 
BSE on the surface P.q=0, which is appropriate for making contact with the mass spectra 
[7]. Its fuller structure is reducible to that of a 3D harmonic oscillator. The ground state 
wave function deducible from this equation has a gaussian structure and is expressible as 
0(g) ~ e~'^^/^^^, where (3 is the inverse range parameter which incorporates the content of BS 
dynamics and is dependent on the input kernel and is given as in [T] [TOl [TT] ) . The structure 
of inverse range parameter /3 in wave function 0(g) is given as [71 [TOl [TT] : 



= (2mim2Mc.yf ) = 1 _ ^ (15) 



We now give the calculation of amplitude and cross section for the process e~ + ^ 
J/ip + rj^ the leading order (LO) of QCD in the next section. We employ the hadron- 
quark vertex functions for 7]^ and J/ijj mesons given in Eq.(ll) and (12) respectively in this 
calculation. 
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III. CALCULATION OF AMPLITUDE AND CROSS SECTION FOR THE PRO- 
CESS e+ + e- ^ J/^ + 

There are four Feynman diagrams in the leading order (LO) of QCD for the process 
e"*" + e~ — > J/^ + f]^. Two of these are depicted in Fig.l. The other two diagrams can be 
obtained by permutations. The details of momentum labeling of the diagram in Fig. la is 
shown in Fig. 2 below. With reference to the momentum labeling in Fig. 2, the adjoint BS 
wave function for rj^ meson can be written down as: 

Qb) = SF{-q2)T''{q,)SF{q,) (16) 

, while for J ftp meson, the adjoint BS wave function can be written as: 

^(Pa, qa) = Sp{-qs)T''{q;,)Sp{q[) (17) 

where qa,qb are the internal momenta of the hadrons J/"^ and rj^ respectively with the 
corresponding hadron-quark vertex functions and F^ given in Eq. (11-12). 




FIG. 1: Two of the lowest order Feynman diagrams for the production of a pair of doubly heavy 
cc mesons in e^e~ annihilation. Other two diagrams can be obtained by permutations. 

Using Feynman rules, one can obtain the amplitude for each of the diagrams in Fig.l. 
The amplitude corresponding to process in Fig. la (described in detail in Fig. 2) is given by 

1 /" 5 

Ml = c5 ^yeeQgl-v{p2)-i ^u{pi) / (fqa(fqbTr[m{Pa,qa)lpSF{qi)ly^{Pb^qb)la]^ (18) 
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FIG. 2: Momentum labeling of the first of the two Feynman diagrams shown in Fig.l 
which can in turn be expressed as, 



^Tr[SF{-q-i)T^{qa)SF{q'i)lpSF{qi)lySF{-q2)^p{qb)SF{qi)la]^ 



(19) 



where c =| is the color factor, the Mendelstam variable s is defined as, s = —{pi +^2)^ 
and eg = 2e/3 is the electric charge of the charmed quark. The momentum relations of the 
quark and anti-quark in the final state are: 

q'l = ^Pa + qa, qs = ^Pa - qa, g4 = ^Pb + qb, q2 = ^Pb - qb, (20) 
and the momenta in the gluon and the quark propagators are given by 

k = q3 + qi = ^{Pa + Pb)-qa + qb, (21) 

q^=q[ + k = Pa + ^Pb + qb, (22) 

respectively. As each of quark momenta in the quark propagators as well as the gluon 
propagator is going to depend upon the internal hadron momenta qa and qb, the calculation 
of amplitude is going to involve integrations over these internal momenta and will be quite 
complex. Hence following |i6J, we simplify the calculation, by employing the heavy quark 
approximation on the quark propagators, where we take the quark masses to be much larger 
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than the internal momenta Qa and qb of the hadrons. In this heavy quark approximation, we 
can use the approximation, << Ma, qb « Mb- Thus we can write, 

Qaib) < Paib) ~ Ma(b), k = ^{P^ + Pb) - qa + qb - \{Pa + Pb) , Ql = Pa + + Qb ^ Pa + \Pb- 

(23) 

With the above approximation, fc^ and q\ are given by A;^ ~ — | and Q'l ~ — | — Tn^c- The 
propagators for the quarks and anti-quarks in momentum space in Eq. (16-17) are given by 
Sp^qi) = _ ^<^^rg7j+m^ ^ whcrc index i labels the quark in the diagram. 

Using Eq.(16) and (17) and the preceding expressions for the gluon and the quark prop- 
agators, the amplitude Mi in Eq.(19) can be written as: 



-2' 2 2-( ^ , , f M ,4 1 r^DiiV.^.(ga)0.(ga)^P^p(?f')0p(9^6) 



—2° f 1 



27(i 2Tii 

(24) 

Here [TR] is the trace over the gamma matrices appearing in the quark propagators in 
Eq.(19). Noting that the 4-dimensional volume element df'q = d^qMda, we then perform 
contour integrations in the complex a-plane by making use of the corresponding pole posi- 
tions [UlIITj. The pole integrations over dq^ = Madda and dq^ = MbdcTb in Eq.(24) can be 
expressed as:;^- r°° ^f"!""" = jJ-^ and tt- = 7=^7^, where values of denominator 

2m J — oo A3A5 iJ{Qa) 2m J—oo A2A4 D{qi,)^ 

functions D{qa) and D{qb) evaluated by contour integration in the complex a- plane are ex- 
pressible as in Eq.(5). After calculating the trace part in the above equation and employing 
the heavy quark approximation on relative momenta given in Eq.(23), one obtains: 



^1 = 32 J ' " ea^^XaSaPbxPaa / d\ld^bN,Np(f),{qa)(l)j,{qb)[v{p2h AP^)] (2^) 

where aem = while = qI/'^'k and is given in Eq.(13). Let's define = / d^4aNy(f)^{qa} 
and = / d^qbNp(j)p{qb), which are values of wave functions at origins of J/ip and i]^ respec- 
tively. 

Thus we can express the amplitude Mi as. 

Ml = ^ ^^"^^'^^ ea^,xaeaPbxPaa[viP2hf,u{pi)]^aU (26) 
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Here Np and A^^, in are the BS normalizers for rj^ and J/ip respectively which are evaluated 
by using the current conservation condition, 

2zP^ = {211)' J d'qTr[^{P,q){^^S^\p,mP,q)S^\-p2)] + (1^2) (27) 

Putting BS wave function ilj{P,q) for a given meson in the above equation, carrying out 
derivatives of inverse quark propagators of constituent quarks with respect to total hadron 
momentum P^, evaluating trace over gamma matrices, following usual steps and multiplying 
both sides of the equation by P^/(— M^) to extract out the normalizer N from the above 
equation, we then express the above expression in terms of integration variables q and a. 
Noting that the 4-dimensional volume element d'^q = d^qMda, we then perform contour 
integration in the complex a-plane by making use of the corresponding pole positions. For 
details of these mathematical steps involved in the calculation of BS normalizers for vector 
and pseudoscalar mesons, see Ref.[Tn]and Ref. [12] respectively, where in the present calcula- 
tion we take only the leading order Dirac covariants i'~f.e and 75 for J/ip and rj^ respectively 
in their respective 4D BS wave functions \1/(P, g). Then numerical integration on variable 
q is performed. The values of BS normalizers thus obtained for J/ip and r] mesons are 

= .0504:GeV-'^ and Np = MlOGeV'^ respectively. 

The total amplitude for the process e'^e'^ — > J /'^rj^ can be obtained by summing over 
the amplitudes of all the four diagrams shown in Fig.l. For that matter the amplitude 
obtained from the first diagram is the same as the amplitude from each of the remaining 
three diagrams in Fig.l. Thus, the total amplitude is 4 times the amplitude from the first 
diagram. The unpolarized total cross section is obtained by summing over various J/^ 
spin-states and averaging over those of the initial state e^e~ . Thus, in the CM frame the 
total cross section, a, is given by 



Ami 


Pf 




327r 


Pi 





j\Y.\^tot?dcose (28) 

spin 



where p/ is the momentum of either of the outgoing particles and pi is the momentum of 
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either of the ingoing particles, which is in turn expressible as, 



4m^ \/s — IQml 

a = — 3 

327r S2 



/" ^ J] I Mtoi |2 dcose (29) 

spin 



where the masses of the leptons are ignored in the calculation. Explicitly | M^f p is given 

by 

1 , , 9^0 4 2 2 6 

\ E I P= ° (-32m; + f + u^)C& (30) 

spin ^ 

where t = — (pi — Pa^ and u = —{pi — PbY are the Mandelstam's variables. Therefore, in 
the CM frame, the total cross section is given by: 



Numerical Results: 

The basic input parameters in the calculation are just four: Co = 0.29, uq = O.lSSGe^, 
QCD length scale, Aqcd = 0.200GeV, and the charmed quark mass, rric = 1.530Gel^ 
flU\ [TT] . The numerical values of inverse range parameter /3 calculated from Eq.(15) are 
= .4989GeV^ and (3^ = ASSSGeV. To calculate the values of (3, for the two hadrons, 
the experimental hadron masses are taken as, Mj/^ = 3.096(^6^^ and M^^ = 2.982Ge^. 
With these parameters the total cross section for the above process at y/s = lO.GGeV is 
calculated to be cr = 21. 75 fb. 



IV. DISCUSSION 

In this paper we have calculated the cross section of the exclusive process of e+e~ — > 
J/'^Tj^ at energy y/s = lO.GGeV in the framework of BSE under CIA [111 HI HI] us- 
ing only the leading order (LO) diagrams in QCD. We find the theoretical value of 
a[e~^e~ — > J/^r^J = 21.75/6, which is broadly in agreement with the Babar's data 
a[e+e~ — > J/'^Vc] = (17.6 ± 2.8 ± 2.1)/6 [T] and the Belle's data, afe+e" — > J/'^Vc] = 
(25.6 ± 2.8 ±3.4)/6 [21 [3]. 

It had been noticed earlier that NRQCD predictions IH [5] for the above process at 
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^/s = lO.GGeV using leading order diagrams alone give cross sections which are much less 
than data [1-3J. Such a large discrepancy between experimental results and theoretical 
predictions has been a challenge to the understanding of charmonium production through 
NRQCD. Many studies were performed to resolve this problem. For instance Braaten and 
Lee [19] first showed that results on cross section are found to improve considerably when 
relativistic corrections are incorporated. Then it was found that to obtain cross sections 
from NRQCD which are consistent with data, one has to incorporate NLO QCD corrections 
[201 EI]- However in these studies it was found that the value of total NLO contribution 
to cross section is nearly twice the LO contribution. In the present calculations under the 
relativistic framework of BSE under CIA, we obtained results for cross sections which are 
in good agreement with data[lH3] using leading order QCD processes alone though we have 
employed the heavy quark approximation [q << M and P ~ M) on the quark and gluon 
propagators on lines of [6j. However we have not made use of the non-covariant heavy quark 
limit as in Eq.(3-4) of Ref.[6]and work with the exact propagators of the quarks constituting 
the two hadrons. This is a validation of the fact that BSE which is firmly rooted in field 
theory and which incorporates relativistic effects within its premises is ideally suited to 
describe not only low energy processes, but even processes at high energies such as high 
energy hadronic scatterings and production processes. 

The approach in this paper is quite different from the approach in Ref.[6] is the sense 
that we employ the framework of BSE under CIA which is a relativistic generalization 
of Instantaneous Approximation (lA) used in the former and has a much wider range of 
applicability as explained in Section II of this paper. Further, to calculate their results, 
Ref . [6] has made use of heavy quark limit on the propagators of all the heavy quarks and 
anti-quarks, where the propagators have been simplified as in their Eq.(3) and (4) which in 
fact is non- relativistic and non-covariant. Also to simplify their calculation of amplitude in 
their Eq.(22), [6] makes use of the heavy quark approximation, in that the propagators of 
quark and gluon are independent of relative momenta qa and qb of the two hadrons since the 
masses of quarks are large compared to their relative momentum. However in our paper, we 
only make use of the above heavy quark approximation {q << M and P ~ M) only in the 
sense of simplifying the integrals involved in Eq.(18)-(19) for amplitude calculation as done 
in [6] , but we do not employ the non- relativistic and non-covariant heavy quark limit on the 
quark and anti quark propagators (as in Eq.(3) and (4) of [6]) and instead work with the full 



14 



quark and anti-quark propagators for the quarks constituting the hadrons. In doing so using 
CIA, we also see that our results on J2s l^iotP in Eq-(30) and cross section in Eq.(31) of our 
paper are not exactly similar to results of [6]. In this regard, we wish to point out that our 
amplitude and cross sectional formulae involve the 4D BS normalizers Ny and Np which are 
calculated in the framework of CIA and whose numerical values are explicitly worked out 
for both the hadrons, J/ip and rj^ as = .0504GeV^~^ and Np = .OAlOGeV'^ respectively. 
These normalizers enter the amplitude and cross sectional formulae through the definitions 
of values of wave functions^^ and C,b of the two hadrons at their origins, which is not so in 
case of Ref. [6]. Further, the input BSE kernel and hence the input parameters employed by 
us and by Ref. [6] are completely different. While we employ "vector" confinement, ie. we 
make use of a common form (7|i^^7|f^) for both one-gluon exchange as well as confining terms 
in the kernel as in Eq.(13) (see[71 [TOjfor details) , Ref.[6] employs a scalar (I'^^^l*^^^) form 
for confinement, while vector form (7[t^^7^^^) for one-gluon exchange. Further the functional 
form of confining potential in [6]is also different from our case. Whereas Ref. [6] uses linear 
confinement (~ r) which is true in case of heavy quark (c,b) systems, we have used a general 
form of confinement potential in Eq.(13) which simulates the effect of linear confinement 
(~ r) for heavy quark sector (large 1711,1712) while retaining harmonic form (~ r^) for light 
quark sector (small mi, 7/12) as explained in section 2. As far as the numerical results on 
cross section in this paper and Ref. [6] are concerned, they are quits close. This may be due 
to the fact that for systems comprising of heavy quarks (c,b), the CIA results may lead to 
lA results when heavy quark approximation is imposed. However this may not be so for 
systems comprising of light quarks (u,d,s). 

However in the present calculation, we used only the first of the leading order(LO) Dirac 
covariants (75 and ij.e)!!! hadron- vertex functions for J/ip and i]^ mesons respectively. They 
were identified as most leading covariants in accordance with our power counting scheme. 
These covariants give maximum contribution to calculation of meson observables such as 
decay constants etc. It can also be seen here that the results on cross section for the process 
e'^e~ — > J/\E' + i]c employing these most leading of the LO covariants brings theoretical 
results close to data [lH3]. We now also intend to see the effect of incorporation of both 
LO and the NLO Dirac covariants (to the vertex functions of these mesons) on cross section 
for the process studied. The contribution from NLO covariants is expected to be much 
lesser than the contribution from LO covariants in line with our recent studies on meson 
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decays [Tot [TTl [15]. And this is more so for heavy mesons comprising of c and b quarks. It is 
expected that the resuhs of cross section will improve further with incorporation of both LO 
and NLO covariants and without employing the heavy quark approximation on the quark 
and gluon propagators. This calculation will be quite rigorous and will be the subject of a 
later communication. 
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